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Two symmetrically coupled populations of N oscillators with inertia m display chaotic solutions 
with broken symmetry similar to experimental observations with mechanical pendula. In particular, 
we report the first evidence of intermittent chaotic chimeras, where one population is synchronized 
and the other jumps erratically between laminar and turbulent phases. These states have finite 
life-times diverging as a power-law with N and m. Lyapunov analyses reveal chaotic properties in 
quantitative agreement with theoretical predictions for globally conpled dissipative systems. 
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Chimera states are remarkable dynamical states 
emerging in populations of coupled identical oscillators, 
where the population splits into two parts: one synchro¬ 
nized and the other composed of incoherently oscillating 
elements [l|. These states have been initially discovered 
in chains of nonlocally coupled oscillators, however they 
can equally emerge in models of globally coupled popula¬ 
tions [ij-lsl. Chimeras have been observed in a repertoire 
of different models [m , and in various experimental set- 
ting s, including mechanical jiaEi, (electro-) chemical 
[l3, llJl and lasing systems [1^, among others. Usually, 
the incoherent oscillators give rise to regular macroscopic 
dynamics which are either stationary, periodic ^o-called 
breathing chimera) or even quasi-periodic [^, 0. Only 
recently, spatio-temporally chaotic chimeras have been 
numerically identified in rings of coupled oscillators [H- 
[l8j| . However, a detailed characterization of the dynami¬ 
cal properties of these states have been reported only for 
phase oscillators with finite-range interactions: in this 
case chimeras are transient, and weakly chaotic [3 [13. 
More specifically, the life-times of these states diverge 
exponentially with the system size, while their dynamics 
becomes regular in the thermodynamic limit. In contrast, 
for globally connected populations, chaotic chimeras have 
so far only been observed in pulse-coupled oscillators [13 
without further analysis. 

In this Rapid Communication, we report the existence 
of various irregular solutions with broken symmetry in 
an experiment with two mechanically coupled popula¬ 
tions of pendula and analyze in depth a model which 
reproduces this complex dynamics. Three pertinent ex¬ 
amples from the experiment, shown in Fig. [T] (a-c), are 
of particular interest. The first two examples represent 
chimeras: in Fig.[T](a), the order parameter of the desyn¬ 


chronized population oscillates quite irregularly, while in 
Fig. [T] (b) it enters a regime of almost periodic oscilla¬ 
tions. Fig. [T] (c) reports a situation where both popula¬ 
tions are irregularly oscillating. We introduce a simple 
model (Eq. ([T])) capable of reproducing all these differ¬ 
ent dynamical behaviors, as it can be appreciated by the 
simulations reported in Figs. [T] (d-f). The model consists 
of two symmetrically globally coupled populations of N 
Kuramoto phase oscillators with inertia. 

The introduction of inertia allows the oscillators to 
synchronize via the adaptation of their own frequen¬ 
cies, in analogy with the mechanism observed in certain 
species of fireflies 2l|. The modihcation of the classical 
Kuramoto model with the addition of an inertial term 
leads to first order synchronization transitions and com¬ 
plex hysteretic phenomena [23 - 13 ]. Furthermore, net¬ 
works of phase coupled oscillators with inertia have re¬ 
cently been employed to investigate self-synchronization 
in power grids 12^31 and in disordered arrays of Joseph- 
son junctions |32l| . In absence of dissipation, Eq. o re¬ 
duces to the Hamiltonian mean-field model - a paradigm 
of long-range interacting systems [13 . 

Our analysis will mainly focus on the solution shown 
in Fig. [1] (a), which we term as an Intermittent Chaotic 
Chimera (ICC). This state exhibits turbulent phases in¬ 
terrupted by laminar regimes, analogous to the one re¬ 
ported in Fig. [T] (b). The third state shown in Fig. [U 
(c) is a Chaotic Two Populations (C2P) state, here the 
erratic dynamics is induced by the evolution of the non¬ 
clustered oscillators belonging to both populations. In 
particular, we show that ICCs are transient states for fi¬ 
nite inertia and system size, whose life-times diverge as 
a power-law with N and m. Furthermore, in the ther¬ 
modynamic limit the intermittent oscillations disappear 
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and the turbulent regime prevails over the laminar one. 
The stability properties of the ICC can be ascribed to the 
universality class of globally coupled systems , which 
are distinct from those reported for chaotic chimeras in 
chains of oscillators M- This result clearly illustrates 
that the stability of chimera states strongly depends on 
the underlying network topology. 
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FIG. 1. (color online) Order parameters (solid black 
cnrve) and (dashed red curve) for the two coupled popu¬ 
lations versus time, (a-c) Experimental measurements; (d-f) 
numerical simulations of the model 0 with m = 10. Ini¬ 
tial conditions are BSCs for (aTidjO) and UCs for (c,f) and 
N = 15. The experiments in panels (a-b) (panel (c)) are car¬ 
ried out with / = 160 beats per minute (bpm) and I = 17cm 
(/ = 184 bpm and I = 25 cm) and experimental time is mea¬ 
sured in seconds. 


Experimental setup. The setup is shown in Fig. [5] and 
it is identical to the one described in 0- The experi¬ 
ment is composed of purely mechanical parts - in partic¬ 
ular, following Q , two populations of non-locally coupled 
metronomes were considered: within each population os¬ 
cillators were coupled strongly, but were coupled more 
weakly to the neighboring population. Metronomes take 


the role of self-sustained oscillators [S^, whose working 


principle is identical to Huygens’ pendulum clocks [36l| . 
except that the escapement in the metronome is driven 
by a spring rather than a mass pulled by gravity. While 
friction inherent to the mechanical elements attenuates 
large-amplitude oscillations toward the unperturbed am¬ 
plitude associated with its unperturbed frequency, small 
pendulum oscillations are amplified by the spring that 
drives the metronome via its escapement mechanism. 
This lends the metronome the characteristics of a self- 
sustained oscillator [^. = 15 identical metronomes 

running with identical frequencies were placed on each 
of two aluminum swings suspended by four rods. The 
strong coupling within one population is mediated by 
the motion of the swing onto which the metronomes 
are attached. As one increases the common frequency 
/ of the metronomes, more momentum is transferred 
to the swing, leading to a stronger coupling among the 


metronomes. A single swing follows a phase transition 
from a disordered to a synchronized state as the cou¬ 
pling within the population increases [s^. The weaker 
coupling between the two swings is facilitated by a pair 
of tunable steel springs, attached to the adjacent rods 
of the opposing swings. The distance of the spring rela¬ 
tive to the pivot can be adjusted: this changes the spring 
lever I and the associated torque, thus effectively tuning 
the spring coupling strength between the two metronome 
populations. 



FIG. 2. (color online) Experimental setup: A sketch of the 
the mechanical system studied iu m, it is composed of two 
swings (A,B) coupled with a spring mechanism, each of which 
is loaded with A = 15 metronomes. 


The motion of swings and metronomes is visualized by 
attaching UV fluorescent spots on swings and metronome 
pendula, and the spot motion is digitally recorded using 
a DSLR photocamera. Subsequently digital data analy¬ 
sis is used to measure the swing and pendulum motion. 

The relative motion of the pendula is reconstructed by 
subtraction of the swing coordinates. Amplitudes and 
phases are then obtained via Hilbert transformation of 
the signal. The phases are used to quantify the level of 
synchronization for each population by using the order 
parameters defined in the following. For exhaustive de¬ 
tails on the experimental setup and methods, see 0 . 

Model and Methods. We consider a network of two 
symmetrically coupled populations of N oscillators. The 
phase of the i-th oscillator in population cr = 1,2 
evolves according to the differential equation 



<t’ = 1 j=l 


( 1 ) 

where the oscillators are assumed to be identical with in¬ 
ertia TO, natural frequency H = 1 and a fixed phase lag 
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7 = TT—0.02. The self- (cross-) coupling among oscillators 
belonging to the same population (to different popula¬ 
tions) is defined as = 0.3 {K^ya' = Ka'a = 0.2), 
with > Ka-a-' as in previous studies on chimera 

states d, [13 • We consider two types of initial condi¬ 
tions: i) broken symmetry conditions (BSCs), realized 
by initializing the first (second) population with identi¬ 
cal (random) phases and frequencies, which may lead to 
the emergence of chimera states, and ii) uniform condi¬ 
tions (UCs) where both populations are initialized with 
random values, and can result in a C2P state 


The collective evolution of each population will be 
characterized in terms of the macroscopic fields (<) = 
i?(°'^(t) exp [j'I'(t)] = The mod¬ 

ulus is an order parameter for the synchronization 
transition being one for synchronous (asyn¬ 

chronous) states. 

The microscopic stability can be measured in terms of 
the associated ordered spectrum of the Lyapunov expo¬ 
nents (LEs) {Xi} i = 1,...,4A^, representing the expo¬ 
nential growth rates of infinitesimal perturbations. The 
dynamics is chaotic whenever the maximal LE Am = Ai 
is positive. For the studied model, presenting a con¬ 
stant viscous dissipative term, the spectrum satisfies the 
following pairing rule [s^l : Xi + A 4 Ar-i+i = —Ijm. 
Therefore the analysis can be limited to the first 2N 
exponents. The Lyapunov spectrum can be numeri¬ 
cally estimated by employing the standard method re¬ 
ported in [ 43 . This amounts to consider for each LE 
Afe, the evolution of a 4A^-dimensional tangent vector 
rW = i = l,...,iV, whose 

dynamics is ruled by the linearization of Eq. (Cl): 


cr' = l 


^ K 


N 


i=i 


( 2 ) 


where = cos {9^^ ^ — 9^'^'^ — 7 ). The orbit and the 

tangent vectors is followed for a time lapse Tg by per¬ 
forming Gram-Schmidt ortho-normalization at fixed time 
intervals At, after discarding an initial transient evolu¬ 
tion Tt- We have employed At = 5 and Tt = 5,000, for 
BSCs we have integrated the system for times 8 x 10"*’ < 
Ts < 3 X 10^ with N = 100,..., 800 and for UCs for 
times 3 X 10"* < Ts < 1 X 10® with N = 100,... ,400. 
The integrations have been performed with a 4‘® order 
Runge-Kutta scheme with time step 5 x 10“^. 

A characterization of the dynamical evolution on short 
time scales can be achieved by considering the probab ility 
distribution function P(A) of the finite time LE A ^|. 
The finite time LE is calculated by estimating the ex¬ 
ponential growth rate of the magnitude of the maximal 
tangent vector 'fA) over finite time windows At, namely 
A = In ||7’*'^^(At)|| where ||7”(^n0)ll = 1- order to 
estimate P(A), we have collected 100,000 data points for 
each system size, obtained from ten different orbits each 


of duration Tg = 100, 000 with At = 10. 

Intermittent Chaotic Chimeras. Starting simulations 
with BSCs at small masses (to < 4), the system is not 
chaotic (as shown in Fig. [3] (a)) and it displays a multi¬ 
tude of coexisting breathing chimeras [3]. Furthermore, 
while the synchronous state RA'i = rA) = \ remains 
stable also in presence of inertia (to > 0), the stationary 
chimeras associated to constant order parameters with 
RA) < rA) = 1 are not any longer observed. For suf¬ 
ficiently large masses, a solution with broken symmetry 
emerges, where one population is fully synchronized with 
RA) = 1, while the other population exhibits wide collec¬ 
tive irregular oscillations in the order parameter RA')(^t) 
between zero and one, as shown in Fig. [1] (d) and Fig. [3] 
(b). These are ICCs and they represent the main subject 
of this Rapid Communication. 



FIG. 3. (color online) (a) Average maximal LE (Ai) vs m for 
N = 100: black circles (red squares) refer to BSCs (UCs). 
(b) Order parameter R (black solid line) and finite time LE 
A (red dashed line) for the chaotic population vs time t for 
N — 200 and m = 10. The corresponding T(A) (solid red 
line) is shown in (c) with the Gaussian fit (black dashed line). 
Inset: probability po versus N. The values (Ai) are obtained 
by following each realization for a time span t = 50, 000 and 
by averaging over 100 different initial conditions. 


As shown in Fig. (b), the erratic oscillations of 
RA) are interrupted by laminar phases, where RA^ stays 
in proximity of one displaying small oscillations. This 
regime is characterized by a large part of the oscillators 
in the chaotic population getting entrained to the syn¬ 
chronous population, apart a few oscillators, which keep 
oscillating with distinct identical frequency, but with in¬ 
coherent phases. An analogous regime is also experimen¬ 
tally observed, as reported in Fig. [1] (b), however due 
to the smaller size of the populations the amplitude of 
the oscillations is larger. Furthermore, by estimating the 
finite time LE A, we show that the laminar phases are 
indeed regular, since they are associated to A ~ 0 (see 
Fig.[3](b)). In particular, the distribution T(A), reported 
in Fig. [3] (c), reveals a clear peak at A = 0, associated to 
the laminar regime, superimposed to a seemingly Gaus- 
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sian distribution. The probability po to observe a laminar 
phase can be estimated by integrating P(A) within a nar¬ 
row interval around A = 0. This probability is reported 
in the inset of Fig. [3] (c) for 10 < At < 1600 and it shows 
a power-law decay with N for sufficiently large system 
sizes, namely N > 50. This is a clear indication that 
the laminar episodes tend to disappear in the thermody¬ 
namic limit. 

To characterize the erratic phase, we give an estimate 
of the average LE A*^*^ restricted to this phase. This es¬ 
timate has been obtained by evaluating the maximum of 
P(A) with a Gaussian fit to the data, once the channels 
around A = 0 were removed. The corresponding data, 
reported in Fig. [4] (a) for various N, reveal a clear de¬ 
cay of A*^*) as l/ln(At). Furthermore, the extrapolated 
value of A^*) ~ 0.022 for N —>■ oo is definitely positive, 
thus indicating that the chaotic state persists in the ther¬ 
modynamic limit for finite m, contrary to what is usu¬ 
ally observed for the Kuramoto model in l^, ^|. This 
logarithmic dependence of the maximal LE A m with the 
system size has been previously reported for globally cou¬ 
pled networks in [3^, where, for dissipative systems, the 
authors have shown analytically that 


Xm{N) = A(=) 


D 

"2 


In(Af) 


O 




(3) 


where is the mean field LE obtained by considering 
an isolated unit of the chaotic population forced by the 
two fields and p^‘^\ D is the diffusion coefficient as¬ 
sociated to the fluctuations of 341. In particular, D 

can be measured by the scaling of the variance of In d(t), 
where d(t) is the modulus of the Lyapunov vector asso¬ 
ciated to the single forced unity. Namely, for sufficiently 
long times it is expected that ([Ind(t) — ~ Dt 

with the average < • > performed over many realizations. 
In the present case, we measured A^'^^ ~ 0.0116(5) and 
D ~ 0.0180(10), thus the expected asymptotic LE should 
be Am(oo) ~ 0.021(1), which is in good agreement with 
the previous reported numerical extrapolation as shown 
in Fig. Idja). This represents the first quantitative verifi¬ 
cation of the prediction ([3]) for a dissipative system with 
continuous time and in particular for an intermittent dy¬ 
namics. 

A more detailed analysis of the stability of this state 
can be achieved by estimating the Lyapunov spectra for 
various system sizes: we observe that the spectrum is 
composed of a positive part made of A^ — 2 exponents and 
a negative part composed of N exponents, see Fig.[4] (b). 
Two exponents are exactly zero: one is always present 
for systems with continuous time, while the second arises 
due to the invariance of Eq. 0 for uniform phase shifts. 
The negative part of the spectrum is composed of an 
isolated LE, quantifying the longitudinal stability of the 
synchronized population, and N — \ identical LEs, which 
measure the transverse stability of the synchronous solu¬ 



FIG. 4. (color online) (a) A^*^ (symbols) versus l/ln(A'') for 
50 < N < 1600 and the corresponding fit, the dot-dashed 
black curve with the shaded area denotes Am(oo) with its 
error bar. (b) Lyapunov spectrum for N = 100, the dashed 
red line is A^'’^ (c) Positive part of the Lyapunov spectra 

for various sizes, the dashed (green) line is A*-'^!. (d) Average 
life-times (r) of the ICC vs N for inertia m = 8 (red squares) 
and m = 10 (black triangles). Inset: (r) vs m for N = 30 for 
ICCs. The blue dashed line in the main panel (inset) refers 
to a power law with exponent 1.60 (2.50). The values of (r) 
are averaged over 200-3000 different realizations of BSCs and 
m = 10 in (a-c). 


tion [i^ . The value of this negative plateau in the spec¬ 
trum coincides with the mean field LE A*^®^ = —0.0266(5) 
calculated for an isolated oscillator of the synchronized 
population, as shown in Fig. S] (b). 


Furthermore, the central part of the positive spectrum 
reveals a tendency to flatten towards the mean field value 
A^®^ associated to the chaotic population for increasing 
system sizes, see Fig. [3] (c) for N = 200,400 and 800, 
while the largest and smallest positive Lyapunov expo¬ 
nents tend to split from the rest of the spectrum. This 
scaling of the Lyapunov spectra has been found to be 
a general property of fully coupled dynamical systems. 
In particular, the authors in 3J, l4^ have shown that in 
the thermodynamic limit the spectrum becomes asymp¬ 
totically flat (thus trivially extensive), but this part is 
sandwiched between subextensive bands containing typ¬ 
ically O(logA^) exponents scaling as in Eq. ([3]) with N. 
We can safely affirm that the chaotic population in the 
ICC reveals properties which are typical of fully coupled 
systems, contrasting with the results reported for chaotic 
chimeras emerging in spatially extended systems 3, 1^ . 


Let us now examine if the ICCs are transient states; 
indeed, we observe for different masses that the chaotic 
chimeras converge to a regular (non-chaotic) state after a 
transient time r. This amounts to the fact that the sys¬ 
tem remains entrapped in a laminar state, which could be 
either a fully synchronous regime or a breathing chimera, 
without returning to the turbule nt p hase. We have mea¬ 
sured the average life times (r) of the ICCs for two 
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masses, namely m = 8 and 10, and various system sizes 
5 < N < 150. These results are displayed in Fig. [4] (d). 
From the figure, it is clear that for TV > 10 one has a 
power-law divergence of the synchronization time with 
an exponent a ~ 1.60(5). The divergence of r is directly 
related to the vanishing of the laminar phases (po —?> 0) 
observable for N —)■ oo. Unfortunately, due to CPU lim¬ 
itations we cannot explore larger system sizes to verify 
that this scaling is present over more decades. However, 
we can safely affirm that these times are not diverging ex¬ 
ponentially with TV as reported in [l6| . This is a further 
indication that our phenomenon has a different nature, 
which is deeply related to the topology presently con¬ 
sidered. Indeed, exponential! 


text of spatially extended systems 


- have been reported for the Hamiltonian version of our 
model [H, Ii3|- Furthermore, we have tested the depen¬ 
dence of (r) on the mass, for one system size, namely 
TV = 30, and we observe that (r) is diverging also as 
a power law of m with an exponent 2.50(5) (see inset 
of Fig. |4] (d)). It is important to remark that regular 
chimeras, appearing for m = 0, are not transient for this 
topology, as we have numerically verified. 

The Chaotic Two Populations State. With UCs, 
the system evolves towards chaotic solutions already at 
smaller masses, namely m > 1, as shown in Fig. [3] (a). 
With these initial conditions the multistability is strongly 
enhanced and many different coexisting states with bro¬ 
ken symmetry are observable, either regular or chaotic. 
By focusing on the chaotic solutions, the so-called C2P 
state reported in Fig. [T] (c),(f) and Fig. [3] (a), we observe 
that in all the cases the oscillators of the two populations 
form a common cluster, characterized by a common av¬ 
erage frequency, plus a certain number of oscillators with 
larger frequencies (see Fig. [5](d)). These states resemble 
imperfect chimeras recently observed in experiments on 
coupled metronomes [Il| and in chains of rotators ^|. 


+ with ||r(’_Ht) = III- The 

time averaged components of the vector are shown in 
Fig.[5](d), from where it is evident that the contribution 
£,i of the oscillators belonging to the common cluster is 
essentially negligible. 
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The C2P states are characterized by a broken sym¬ 
metry since the dynamics of the two populations takes 
place on different macroscopic chaotic attractors, as is 
clearly observable in Fig. [5] (a) and (b). In particu¬ 
lar, the differences in the oscillation amplitudes of 
and are due to the different number of oscillators 
contributin g to the common cluster in the two popu¬ 


lations [2^ l23|, l27|. The multistability is also reflected 


in the associated Lyapunov spectra; three examples are 
shown in Fig. [5] (c), their shapes are extremely differ¬ 
ent presenting different numbers of positive and negative 
LEs. However, a general result is that the oscillators 
contributing to the chaotic dynamics are the ones out of 
the common cluster. As shown in 0 , the contribution 
of each oscillator i to the chaotic dynamics can be mea¬ 
sured in terms of the corresponding squared component 
of the normalized maximal Lyapunov vector namely 


i/N 


FIG. 5. (color online) (a) Order parameters (black solid 
line) and (red dashed line) vs time; (b) imaginary vs real 
part of the complex fields and (c) Lyapunov spec¬ 

tra obtained for 3 different realizations of UCs. (d) Averaged 
components U (red diamonds) together with the correspond¬ 
ing averaged frequencies of the oscillators chi = ^ (black 
circles) for both populations. All the data refer to C2P states 
and to m = 9, the system sizes are N = 200 in (c) and A = 50 
in (a),(b),(d). The time averages have been performed a time 
T ~ 5 X 10'‘. 

Conclusions. We have shown that a simple model of 
coupled oscillators with inertia can reproduce the erratic 
behaviors observed in our experiment of two coupled pop¬ 
ulations of mechanical pendula. The presence of iner¬ 
tia is a distinctive ingredient to observe the emergence 
of chaotic regimes, like ICCs and C2Ps. The detailed 
characterization of the ICC dynamics reveals that its 
chaotic properties can be interpreted in the framework 
of fully coupled dissipative systems 


34| . However, our 


study extends the validity of the results reported in [3^ 
to networks with inhomogeneous coupling displaying in¬ 
termittent dynamics. Together with the results reported 
in for a ring geometry, this clearly indicates that 

the stability properties of chaotic chimeras strongly de¬ 
pend on the underlying network topology. It would be 
extremely challenging to investigate if the topology also 
influences the stability of non-chaotic chimeras. 

Our dissipative model is quite remarkable, since it dif¬ 
fers from a Hamiltonian model only by a constant dis¬ 
sipative term proportional to I/m that vanishes in the 
limit of large inertia. This suggests that for sufficiently 
large m dynamical properties characteristic of conserva¬ 
tive models should be observable. Indeed, the Lyapunov 
spectrum exhibits a pairing rule [s^ similar to that of 
Hamiltonian models. Furthermore, ICCs are metastable 
states, whose life-time diverges as ~ TV^ ®, in analogy 
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to quasi-stationary states in the Hamiltonian mean-field 
model [ 3 ^. 
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